clc; clear;
cd 'data'

lambda_data = csvread('lambda.csv');
rho_data = csvread('rho.csv');
Y_data = csvread('Y.csv');

cd ..

I = 44;

alpha = 1; 
theta = 2.5; sigma = 3;
YEAR = 15;

lambda = reshape(lambda_data(:,YEAR),I,I);    
rho = 1-reshape(rho_data(:,YEAR),I,I);   
rho_autarky =  max(0,1-(1-diag(rho))./diag(lambda));
 
Y = Y_data(:,YEAR);

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = (Gini_autarky./Gini);

%---------------------------------%

theta = 5; sigma = 3;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

%---------------------------------%

theta = 8; sigma = 3;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

%---------------------------------%

theta = 5; sigma = 4;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

%---------------------------------%

theta = 8; sigma = 4;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

%---------------------------------%

theta = 5; sigma = 5;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

%---------------------------------%

theta = 8; sigma = 5;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

%---------------------------------%

theta = 8; sigma = 6;

Gini = alpha.*(sigma-1)./(sigma*theta) - 2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1))...
    .* sum(lambda.*(1-rho).*repmat(Y',I,1)./repmat(Y,1,I),2);
Gini_autarky = alpha.*(sigma-1)./(sigma*theta) -  2.*alpha.*(theta-(sigma-1))./(sigma*theta).*(sigma-1)./(4*theta-2*(sigma-1)).* (1- rho_autarky);
 
Gini_change = [Gini_change, (Gini_autarky./Gini)];

Table = Gini_change;

clearvars -except Table


